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Abstract: The primary goal of this paper is to characterize solutions to coupled 
reaction-diffusion systems. Indeed, we use operators theory to show that under 
suitable assumptions, then the system given by 

u t = M A u + F(u) 

have solutions. As applications, we consider a mathematical model arising in 
Biology and in Chemistry. 



1 Introduction 

The primary goal of this paper is to make an investigation on a particular 
type of partial differential systems, that is, the strongly coupled reaction- 
diffusion system. Recall that early investigations on this problem are due 
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to Fichera, see, e.g., [8], and recently by Amann, see, e.g., [1, 2, 3]. Also, 
recall that in most of publications on the proposed problem, the diffusion 
matrix M is supposed to be diagonal and sometimes with positive entries. 
Obviously such an assumption cannot be applied to some interesting cases 
arising in several fields such as in Biology, Chemistry, and Ecology. In this 
paper, we consider the general case, that is, we assume that M is any matrix 
without any restrictions on its entries. However it will be shown that if all 
eigenvalues of M belong to S = {z e C : 3t,ez > 0}, then the reaction- 
diffusion systems have a unique solution. The main idea of our investigation 
is based on operators theory, especially unbounded normal operators and 
related semi-groups of contraction. The strongly coupled reaction-diffusion 
system is defined as 

u t = MAu + F(u), (1) 

where M is a d x d real matrix and F : H d i— > R d is of class C 2 . Throughout 
this paper we will assume that d is even; the case where d is odd will be 
investigated elsewhere. However, the author expects to use a similar method 
as in this paper. 

As stated above, we propose to solve Eq.(l) using unbounded normal 
operator method. Under appropriate hypotheses, we prove that Eq.(l) has 
a unique solution. 

Let us consider the needed background for it. Let A (not necessarily 
bounded) be a normal operator in the (complex) space Hilbert H ( A is 
densely defined, and closed such that A A* = A* A). Using the spectral 
theorem for unbounded normal operators, see, e.g., [4, or II, pp. 348-355], 
it is well-known that A = A\ — iA 2 , where A\ and A 2 , are respectively 
the real part and minus the imaginary part of A. Notice that A^s are 
respectively self-adjoint operators. It is also well-known that if we assume 
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(So) 



that AkS are nonnegative self-adjoint operators, then iA is an m-accretive 
operator. Thus (iA) is the infinitesimal generator of a contraction semi- 
group, see, e.g., [4, or 10 Corollary 4.4, p. 15]. To examine Eq.(l), we study 
the linear part of it, that it, the diffusion operator L = MA in the Hilbert 
space H = [L 2 (Vi)} d where Vl is a bounded subset of R p with a smooth 
boundary. In fact, we need to study the following problem 

u t = -MAm 

. u \an = 

More generally, let A, B, C, and E be unbounded normal operators in 
H, and consider the following system 

Ut = Au + Bv 
v t = Cu + Ev 
u(0, x) = Uo(x) 
v(0, x) = vo(x) 

Such a system is equivalent to the following: 

X t = T X with X = (u, v) 
X(0,x) = X (x) = (u (x),v (x)) 
where the unbounded matrix operator T is defined by 

T : [D(A) n D{C)\ © [D{B) n D(E)} ^H ®H 

and 



(Si 



(S2) 




T = 



Set T = —iT. Thus, the unbounded operator T is defined by 

/ 

f = 

iC %E 



i A iB 



\ 



It follows that (Si) <^> (S 2 ) O (S 3 ), where (S 3 ) is given by 
(S3) < 



X t = -iT X with X = (u, v) 
X(0,x) = Xq(x) = (u (x),v (x)) 



The next step is to show that T is a normal operator, which implies that 
— T is an m-accretive operator. 

2 Diffusion Equation 

In this section, we show that the general problem (Si) admits a unique 
solution under appropriate hypotheses on A, B, C, and E. As particular 
case, the equation (S ) will be considered. 

Theorem 2.1.- Let A, B, C , and E be unbounded normal operator in 
the Hilbert space H . Assume the following conditions hold true 

1. N(A) = {0}, and N(E) = {0} 

2. D(A) n D(C) and D(B) n D(E) are dense in H 

3. A~ X B and C~ X E are closed operators in H 

4- A, B, C , and E commute each other 

where N(A) and N(E) are respectively the Kernels of the operators A and 
E. Then the matrix operator T is normal in H © H . 

Proof. Recall that the matrix operator T is defined in H®H by, D(T) = 
[D(A) n D(B)} © [D(C) n D(E)\, f(u, v) =< i(Au + Bv),i(Cu + Ev) > 
y(u,v) G D(T). Therefore, it is a densely defined operator in H © H, 
according to assumption (2). Let (u n ,v n ) be a sequence in D(T) such that 



(u n , v n ) converges to (u, v) and T(u n , v n ) converges to irj) in H © H. In 
other words, Au n + £?i> n and Cii n + Ev n converge to £ and rj respectively. 
Since the kernel N(A) = {0} (according to (1) ), then A~ x Bv n — > A~ x ^—u. 
Now, since A~ l B is closed then v e D(A- l B) = D(B) and A~ l Bv = 
- u. In addition u e D(A), since u = A' 1 ^ - A' x Bv e Thus, 
it easily follows that (u, v) e Z)(A)©D(B) and A-u+_Bv = ^. Using a similar 
argument yields (u, v) G -D(C) © and Cm + Ev = r/. Therefore, T is 

a closed operator. Now, we have 

/ 



A*A + C*B B* A + E*B 
A*C + C*E B*C + E*E 



and, 

\ 



) 



AA* + CB* BA* + EB* 
AC* + CE* BC* + EE* 
Clearly AA* = A* A, and EE* = E*E( A and E are normal operators). 
Let show that C*B = B*C. A similar argument can be used to show that 
BA* = AB*, E*B = EB*, and C*E = CE*. Let us write C = C x -iC 2 and 
B = B x — B 2 as stated in the introduction of this paper. Thus C* = Ci+iC 2 
and B* = B 1 + iB 2 . Therefore C*B = (d^ + C 2 B 2 ) + i{C 2 B x - C X B 2 ) 
; in the same way B*C = (B-yCi + B 2 C 2 ) + i(B 2 Ci - BiC 2 ). Now, since 
BC = CB (according to (4)), we have C p B q = B q C p for p,q — 1,2. It 
follows that C*B = B*C. In summary, we have f*T = ff*. The proof is 
complete. 

Corollary 2.2 Under previous assumptions. The operator T can be 
decomposed as T = 7\ — iT 2 where T\ and T 2 are respectively self-adjoint 
operators. Assume that both 7\, T 2 are nonnegative operators. Then the 
operator —T = %T is m-accretive. In addition 5i admits a unique solution. 



Prof. Since T is a normal operator and that both its real and minus 
imaginary parts are nonnegative self-adjoint operators, then iT = —T is m- 
accretive, see, eg., [10, Corollary 4.4, p. 15]. It follows that the system (Si) 
admits a unique solution. Equivalently both (S 2 ) and (S 2 ) admit unique 
solutions. 

We apply previous results to the problem (S ). Assume that d = 2n and 

set 

A = —Mi A, B = -M 2 A (2) 
C = -M 3 A, and E = -M 4 A 
where M k (k = 1,2,3,4) is anx n-matrix. Consider the following problem 



(S. 



4 



X t = -MA X with X = (u, v) 
X(0,x) = X (x) = (u (x),v (x)) 

where MA : [H 2 (tt) n ^o(^)] 2n ^ [^ 2 (^)] 2n is defined as 



MA 



^ Mi A M 2 A 



V 



M 3 A M 4 A 

Set T = -MA and — T = iT. Therefore T is defined as 

( —iM 1 A -iM 2 A 



T = 



v -iM 3 A -iM 4 A 



According to theorem 2.1, T is a normal operator. We have the following 
result. 

Theorem 2.3 Under previous assumptions. Assume that M is not the 
zero matrix. In addition if all eigenvalues of M belong to the set S = {z G 
C : 3t,ez > 0}. Then the problem (S 4 ) admits a unique solution. 
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Proof. The operators A, B, C, and E given in Eq.(2) satisfy assumptions 
(l)-(2)-(3)-(4) of the theorem 2.2. It turns out that T is a normal operator. 
Now let us show that if all eigenvalues of M belong to S = {z e C : ^tez > 
0}, then iT = MA is m-accretive. Let Ai, Ai, A r ( ki + k 2 + ... + k r — 2n) 
be eigenvalues of M. Following the Jordan decomposition method for M, 
it is well-known that MA can be decomposed as 



ma = n 



J fcl (Ai)A 

J fe2 (A 2 )A ... 



\ 
where II is a nonsingular matrix, and 








n- 1 



J fer (A r )A J 



( 



J fe (A)A 



AA A ... 
AA A .. 



k > 2; Jx(A) = [AA], k = 1 



y ... AA j 

Thus, all eigenvalues of M belong to S = {z e C : $tez > 0} if and only 
if iT = MA is m-accretive. Therefore, if all eigenvalues of M belong to 
S then iT = MA generates a contraction semi-group. In such a case (S4) 
admits a unique solution. Equivalently (S ) admits a unique solution, since 
(Sq) is a particular case of (S4). 



3 Coupling Problem 

Assume d — 2n and consider the coupling of the diffusion with the reaction 
term, that is, the system given by Eq.(l). We define the following operators 



D(T) = [H^Q)nH\Q)] d } 
Tu = MAu VueD(T) 

where M is the d x rf-matrix given in the Eq.(l). 

D(R) = {ue [L 2 (Q)] d : F(u) G [L 2 (Q)] d } 
Ru(t,x) = F(u(t,x)) a.e u G D(R) 

We will make the following hypotheses 

(H ) all eigenvalues of the matrix M belong to S = {z G C : ?R,ez > 0} 

(Hi) Assume the operator R is m-accretive, and that G i?(0) 
We have the following. 

Theorem 3.1. Under assumptions (H ) and (H ). Then the Eq.(l) 
admits a unique solution. 

Proof. The main idea is to show that the nonlinear operator given by 
T + R is m-accretive in H = [L 2 (Q)] d . Consider Yosida's approximation for 
R. It is defined as 

Rx = j[I-(I + XR)- 1 ], A > (3) 

It is well-known that R\ is m-accretive and that —-Lipschitz in H. Now 

A 

consider the following equation 

eu\ + Tu\ + R\U\ = v, and u\ = 0, on dVl (4) 

Since T + R\ is m-accretive, see, e.g., [5], then Eq.(4) admits a unique 
solution u\ G D(T) for any v G H, and A > 0. We also know the family 
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1 



( u x)x>o is bounded by - \\v\\ H . Using the fact R\ is m-accretive, R\0 = 0, 

e 

and by integration by parts it easily follows that 



Thus, multiplying Eq.(4) by Tu\, and from Eq.(5), it turns out that (Tu\), 
and (R\U\) are bounded. From the compactness embedding, [H 2 (Q)] d 
[L 2 (n)} d , and the fact (u x ), (Tu x ), and (R x u x ) are bounded, it turns out 
that: u\ strongly converges to u, (Tu\) weakly converges to £, and (R\U\) 
weakly converges to 77, as A approaches to in H. Since T is closed, then 
Tu = £. Since R is m-accretive, then Ru = 77, see, e.g., [5]. In summary 
T + R is m-accretive under assumptions (Hq) and (Ho). Therefore the 
algebraic sum(see, e.g., [6]) (T + R) generates a nonlinear contraction semi- 
group, that is the Eq.(l) admits a unique solution. 



In this section, we consider a model considered in [9]. The problem we will 
study represents a mathematical model describing various chemical and 
biological phenomena. In [9], Lyapunov functionals have used to prove a 
global existence of unique solutions. Here, we use the method described 
above to prove that the given problem admits a unique solution, under 
suitable assumptions. 

Our model is described as 




d 



(5) 



4 Applications 



(M) 



Ut — aAu — j3Av = —af(u, v) in (0, 00) x Q 
v t — 7 Au — aAv = pf(u, v) in (0, 00) x Q 
u u = v u = on (0, 00) x dQ 
u(0,x) — u (x), v(0,x) — v o (x) in Q 
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where u(t) and v(t) represent either chemical concentrations or biological 
population densities, Q is a bounded open subset of class C l in R n , u v (re- 
spectively v v ) denotes the outward normal derivative on dfl, and a, [3, 7, p, 
and a are positive constants . In [9], the following hypothesis is made 



2a > (P + 7) 



(6) 



The (M) can be expressed as 



(AO 



(u t ,v t ) = MA(u,v) +F(u,v) 
(uu,v v ) = (0,0) 

(u(0,x),v(0,x)) = (u (x),v (x)) 



where F(u,v) = (—af(u,v) , pf(u,v)), and 



M = 



a (3 
7 a 



It is obvious to see that all eigenvalues of of the diffusion matrix M are 
given as, EV(M) = {a + v^7 , ot — \JJh}. Since all eigenvalues of M are 
nonnegative, then 

« > V^7 (7) 

Clearly Eq.(6) implies Eq.(7). Indeed, ^((3 + 7) > yfjTy. Therefore, 
instead of considering Eq.(7), we will only assume that Eq.(6) holds. 
Consider the Hilbert space H = [L 2 (Q)] 2 and set 



D(T) = [H 1 (Q)nH 2 (Q)f 



and 



T(u,v) = M(Au,Av) 
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In the same way, define 



D(R) = {(«, v) G [L\Q)] 2 : (-*/(«, v),pf(u, v)) G [L 2 (0)] 2 } 



and 



R(u,v) = (—af(u(t,x),v(t,x)) , pf(u(t,x),v(t,x)) a.e u,v & -D(-R) 
We will make the following assumption 



For instance from the fact that R is accretive, the following holds 



More generally, assume that / is given such that R is a nonlinear m-accretive 
operator in [L 2 (f2)] 2 . Thus, we have the following. 

Proposition 4.1. Under Eq.(6), and Eq.(8), then the problem de- 
scribed in (M) admits a unique solution. 

Proof. Obvious as consequences of Theorem 2.3 and Theorem 3.1. 
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